High energy neutrinos from charm in astrophysical sources 
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Charm production gives rise to a flux of very high energy neutrinos from astrophysical sources with 
jets driven by central engines, such as gamma ray bursts or supernovae with jets. The neutrino flux 
from semi-leptonic decays of charmed mesons is subject to much less hadronic and radiative cooling 
than the conventional flux from pion and kaon decays and therefore has a dominant contribution at 
higher energies, of relevance to future ultrahigh energy neutrino experiments. 
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I. INTRODUCTION 



Large underground or underwater experiments like Ice- 
Cube and KM3NeT [| are designed with the goal of 
observing high energy neutrinos produced in astrophysi- 
cal sources. The highest energy neutrinos, with energies 
of 10 9 GeV and higher may be observed in radio detection 
experiments Q, and with an even higher energy thresh- 
old of 10 12 GeV with acoustic detection experiments 
We consider astrophysical sources driven by a relativistic 
jet outflow, accelerated by a central engine such as a black 
hole 0, H[ . Shock accelerated protons in the jet outflow 
may give rise to a high-energy neutrino flux 0. These 
neutrinos are potentially produced in hadronic interac- 
tions: proton-proton interactions produce charged pions 
and kaons which subsequently decay into muons and neu- 
trinos. Above the threshold for A + production, proton 
interactions with ambient photons also produce charged 
pions, and at higher energies, kaons. The relative im- 
portance of the pp and pry contributions to the neutrino 
fluxes depends on the characteristics of the astrophysical 
environment. Several types of astr ophysic al sources have 
been studied in e.g. 0, H i EE El El El 0, EE El S3 • 

High energy pions and kaons are relatively long-lived 
and therefore subject to both hadronic and radiative 
cooling before they decay, which downgrades the neu- 
trino energies. Charm production and decay in astro- 
physical jets is also a source of neutrinos In this 
paper, we show that production of charmed mesons in 
pp collisions gives a large contribution to the neutrino 
flux at the highest energies, since high energy charmed 
hadrons {D^, D°) have short lifetimes and therefore pre- 
dominantly decay before they interact. Moreover, since 
the amount of radiative cooling scales as m -4 , the larger 
masses of the charmed hadrons lead to less cooling. The 
neutrino flux from charm is therefore less suppressed up 
to higher energies. Even though the production cross 
section is orders of magnitude smaller than for pions and 
kaons, neutrinos from charm decays become the domi- 
nant contribution at high energies. The energy at which 



charm begins to dominate depends on the detailed prop- 
erties of the astrophysical source, and the extent to which 
charm contributes depends on the maximum energy of 
the accelerated protons in the jet. 

An example of an astrophysical source model is the 
slow-jet supernova (SJS) model, proposed by Razzaque, 
Meszaros and Waxman (RMW) 9]. It is characterized 
by a mildly relativistic jet propagating in a collapsing 
star. This jet does not emerge from the source, and is 
sometimes referred to as a 'choked jet.' This environment 
has a large optical depth [ljl EH , so neutrinos may be 
the only high energy signals to emerge. On the other 
hand, the jet in a Gamma Ray Burst (GRB) (see e.g. @ 
for a review) is highly relativistic and the optical depth 
might be such that photons escape during the burst and 
do not thermalize. Other GRBs may be choked, so that 
only neutrinos escape 

We present here a new analysis of the neutrino flux for 
two types of source environments: the SJS and a GRB 
model, both of which have large magnetic fields in the 
jets. Our treatment of the proton-proton collisions ac- 
counts for charmed meson production, as well as pion 
and kaon production. The shorter lifetimes of charmed 
mesons allow this channel to potentially dominate the 
neutrino flux from pp collisions. 



II. ASTROPHYSICAL SOURCES 

Cooling times for protons and mesons depend on the 
proton and photon densities in the jet, which are deter- 
mined from the characteristics of the environment. The 
internal shock radius of the jet is given by Tj — 2TjCt v , 
where Tj is the bulk Lorentz gamma factor of the jet 
outflow and t v is the variability time scale of the central 
engine. The proton and electron number densities can be 
calculated in terms of the jet luminosity, Lj, and the jet 



opening angle, 9j, as 



L : 



2^ 2 r 2 r 2 m pC 3 
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The primes denote quantities in the frame comoving with 
the jet. The average magnetic field generated in the jet 
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is given by 



B' 



1/2 



(2) 



where es is the fraction of kinetic jet energy converted 
into magnetic fields. The comoving energy density of 
photons in the jet is 



(3) 



where e e is the fraction of the energy radiated into pho- 
tons. 

Depending on the optical depth, in some astrophysical 
sources photons may be thermalized. For a thermal pho- 
ton distribution with temperature T", the photon number 
density is given by n' oc (fcT') 3 , leading to 
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In the astrophysical environments where photons are not 
thermalized, thephoton energies follow a broken power- 
law spectrum [18] . characterized by the break energy 
E' l, taken as a free parameter, and U' 



A. Slow Jet Supernova Model 



The RMW model of this system has a bulk jet Lorentz 
factor Tj = 3. The characteristic time variability of 
the source is t v = 0.1 s, with a jet luminosity of Lj = 
3 x 10 50 erg/s and jet angle of 9j ~ 1/3. We take the 
energy fractions e ei s =0.1. With these parameters, the 
estimated size of the jet is on the order of 10 11 cm, below 
a typical stellar radius for a pre-supernova star. Due to 
the large optical depth, photons thermalize. 



B. Gamma Ray Bursts 

We also consider a generic case of a GRB. There is a 
wide range of possible parameters and types of progeni- 
tors of GRBs, and the choice of parameters we make here 
is to be considered as just one example of a GRB-like jet. 
The jet gamma factor can vary over a large range of val- 
ues, and the medium can be optically thick or thin. We 
choose the case of a GRB where the photons are not in 
thermal equilibrium and we take the values Tj = 100, 
Lj ; = 3 x 10 50 erg/s, 6j = 0.1, t v = 10" 3 s. We fix 
the photon break energy to be E' b = 2.5 keV [l9[ and 
the spectral indices to be —1 below and —2.25 above the 
break energy. We take the same energy fractions, e e ,B, 
as for the SJS model. 

In Table [T] we show the number densities of protons 
and photons, the values of the magnetic field and the 
characteristic photon energy, in the comoving frame, for 
the SJS and GRB models. 



TABLE I: The jet bulk Lorentz factor Tj, and the comoving 
number densities of protons n' p and photons ri 1 , average pho- 
ton energy EL,, and magnetic field in the jet B' for the slow-jet 
supernova (SJS) and gamma ray burst (GRB) models. 



Source 


Tj 


n' p [cm 3 ] 


B' [G] 


EL, [keV] 


nL, [cm -3 ] 


SJS 


3 


3.6 x 10 20 


1.2 x 10 9 


4.5 


2.8 x 10 24 


GRB 


100 


3 x 10 16 


1.1 x 10 7 


2.5 


1.1 x 10 21 



III. TIME AND DISTANCE SCALES 

The neutrino flux is determined by the flux of acceler- 
ated protons, the production of mesons, meson interac- 
tions and meson decay. This all occurs in an environment 
of magnetic fields, photons, and a non-relativistic baryon 
density. To begin with, we consider the characteristic 
time scales and distance scales relevant to the acceler- 
ated protons. Since many of the processes are relevant 
to mesons as well, we also indicate the application to 
mesons. 

Protons are accelerated by Fermi shock acceleration. 
The acceleration time is 



3 x 10- 12 s K °-3 f E 'p\ 
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B' 
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where the the parameter re is inversely related to the dif- 
fusion coefficient. The details of the shock acceleration, 
in particular the orientation of the magnetic field relative 
to the shock, determine the value of diffusion coefficient 
[20| . Using the standard choice of es =0.1 and re = 10, 
the acceleration time converted to a distance is 



t acc 



10" 



(3>\ 



B 1 
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(6) 



A value of re = 10 taken here could be reduced by a fac- 
tor of 10 or more depending on the orientation of the 
magnetic field (20l . [2l| . We comment below on the impli- 
cations of a smaller value of re. 

Until proton cooling times (or proton interaction 
lengths) are smaller than the acceleration time, the pro- 
ton energy spectrum within the jet is characterized by a 
power law which we take to be 4>n ~ (E' p )~ 2 . We neglect 
shock acceleration of the charged mesons to simplify the 
evaluation [TTj |. 

The time scales for proton cooling depend on the 
comoving energies, particle densities, cross sections, 
and magnetic fields. Generically, the cooling time is 
Col = E'/\dE'/dt'\. One can approximate \dE'/dt'\ ~ 
cn'aAE' , where AE' is the particle energy loss and a is 
the cross section of the hadron scattering with particles 
with comoving number density n' . For scattering on pho- 
tons, we use (n'av), averaged over the photon spectrum. 

There are several processes that contribute to the pro- 
ton and/or meson cooling. The effect of the energy loss, 
accounting for the weighting by the incident flux, can be 
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described by a cascade equation written in terms of Z- 
moments [221 ] - For the proton flux (<^jv) the propagation 
over distance x' in the co-moving jet frame is given by 



dtp 



TV 



dx' 



4>N 



cool 



<PN „had 

flhad ~*~ NN ^had 
"N 4 TV 



(7) 



where 



((TppTip) -1 describes interactions with tar- 
get protons. The ellipsis signifies other cooling mecha- 
nisms, such as p7 reactions, and electromagnetic cooling, 
to be discussed below. The Z-moment is defined by 



dXEXg' 



_l driN- 



• M 



dxi 



(8) 



where the flux of N is <j> N ~ (£")~ Q , x E = E' M /E' N 
and dn/dxE describes the cross section normalized en- 
ergy distribution of the meson M produced by N (or 
from TV decay). 

Using the comoving proton number density from Ta- 
ble [J and the cross section a pp ~ 5 x 1CP 26 cm 2 , the 

5.5 x 10 4 cm for the SJS 
6.8 x 10 s cm for the GRB. 
We used these same interaction lengths for the hadronic 
interaction of the mesons. Using the exact value for the 
inelastic cross section for mesons does not affect the re- 
sult because the hadronic cooling time for mesons is much 
larger than its radiative cooling times. 

The distance scale corresponding to the hadronic cool- 
ing time is 



interaction length is £^ d 
model of RMW, and 



j had 



/'had 

-/V 
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yhad 



.had 
CZ N . 



(9) 



The Z-moment takes the inelasticity into account, and is 
equivalent to setting AE' /E' ~ 1 — Z^f^ for this hadronic 
scattering process. We will use Z^^ — 0.5 [23j ]. 

In addition to hadronic interactions, protons can lose 
energy via photon interactions such as inverse Compton 
(IC) scattering, where the proton loses energy to the am- 
bient photons, or through inelastic pj interactions if the 
proton energy is above the A + threshold. For the IC pro- 
cess, as long as the center of mass energy is small com- 
pared to the proton mass, one can use the Thomson cross 
section, <tt p = <Tt (to 2 / 'm 2 ) . This approximation is good 
for protons scattering off kcV photons for energies up to 
E' p ~ 10 5 GeV. Above this energy, IC becomes irrelevant. 
The effective scattering length for inverse Compton scat- 
tering is 



TIC _ 

n N — 



o~~4 4 

Am p c 



4a T m 2 e E>U> 



(10) 



This effective scattering length is rescaled by (tum /rn p ) A 
for mesons. 

For E^ = 5 keV, the threshold energy for A + produc- 
tion in p7 interactions is E' p th = 2 x 10 5 GeV. For pj 
scattering, the averaged reaction rate is given by, 




FIG. 1: Proton cooling times for SJS and for GRB in the 
frame co-moving with the jet. 



(n av) — 



dE'-r E a J ds{s~m 2 p )a pl {s), 

(11) 

where /3' p is the proton beta factor, s is the center-of-mass 
energy squared, and fU/(E')dE' is the photon number 
density in the energy range between E^ and EL + dE^ . 
The cross section cr p7 is given by the resonance plus con- 
tinuum multiparticle production contributions given in 
Ref. [13]. The photon distribution rL is thermal for the 
SJS and a power law for the GRB j7|. In determining 
the cooling time, t', we take AE' /E' w 0.2 1 . 

Finally, there is a cooling distance associated with syn- 
chrotron radiation because of the magnetic field in the 
jet. This distance scale is 



L 



N 



a T m 2 E'B r - 



(12) 



Again, meson interaction lengths can be obtained by scal- 
ing by (m M /m p ) 4 . 

The cooling times and the acceleration time for protons 
are shown in Fig. [TJ For the proton flux, until cooling be- 
comes important, acceleration dominates. For both the 
SJS and GRB example models, cooling through proton- 
photon interactions eventually dominates the accelera- 
tion. The crossover between acceleration time and t pl 
in Fig. 1 occurs in an energy regime where t pl is nearly 



The Bethe-Heitler process p7 — > pe + e is less effective than 
inelastic p-y above the threshold [9| and we will therefore neglect 
it. 



4 



constant, although there is strong energy dependence at 
lower energies. 

When cooling dominates acceleration, the proton flux 
is cut off. Details of the cutoff depend on the energy 

tendence and mean inelasticity of the interaction [2ll . 
Generically, we can write 

<j> N (E) ee <p N (E)f N (E) = AE- 2 f N (E) . (13) 

We apply a smooth cutoff Jn{E) of the form of Eq. (20) 
in Ref. [2l[ with a spectral index T = 2 and the parameter 
6 = 1: 



Jn{E) = 1 



E 



E n 



exp 



(T-l) 



E 



E n 



The energy at which the cross-over of acceleration and 
cooling times occurs, in the co-moving frame, is E max ~ 
5.2 x 10 5 GeV for the SJS and £ max ~ 3.6 x 10 6 GeV 
for the GRB model. Our result for the SJS model is 
in agreement with previous estimates of proton cooling 
times dlH. 

Fig. [2] shows the meson cooling times and the decay 
times in the frame co-moving with the jet for both the 
SJS and GRB models. For mesons, the decay length and 
decay time includes the gamma factor (in the comoving 
frame), e.g., L^ c = (E' M /itimc 2 )ctm ■ The proper decay 
lengths ctm are listed in Table HT1 Pion and kaon cooling 
times are in agreement with previous estimates [H, [To| . 
One can see from Fig. [5] that for D mesons, the cross 
over of the decay time with the cooling time (hadronic or 
radiative) occurs at much higher energies than for pions 
and kaons. 



IV. MESON AND NEUTRINO FLUXES 

Neutrinos are produced through meson production fol- 
lowed by meson decay. It is convenient to define the com- 
bined cooling length Lm for meson M through 



N 



[(L^)" 1 + (Z^)" 1 ] 



(14) 



M 



including the mass rescaling factor. An effective length 
for a given meson which also accounts for the decay is 
(iff)" 1 = (Lm)- 1 + {L% c )~ l . With this definition, the 
equation describing the propagation of mesons M in the 
jet is 



d<t>M 
dx' 



reft 

M 



Znm 



4>N 

/>had 
N 



7 7 <t>N 
^NM »7 



(15) 



where Znm are the Z-moments for pp meson production 

J NM 



and Z%ri„ are the Z-moments for pj meson production. 



We have here assumed a constant density that does not 
vary with x' . Finally, the neutrino flux from meson M is 
obtained from the meson flux by integrating the equation 



dx' 



<Pm 

t dec 
M 



(16) 
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FIG. 2: Meson cooling times for SJS and GRB in the frame 
co-moving with the jet. 



where Zmv is the decay Z-moment listed in Table [TTJ In 
the limit x' — * oo we have 



ME 1 ) = z 



77 «had 
NM N 



T doc 
M 



/?had 
N 



■ (-If 



<Pn(E') . (17) 



This has the correct limiting behavior, e.g., when decays 
dominate the meson effective length and the energy 
is below the pr) cross section threshold. In this limit, 

ME') = Z Mv Z NM <t>N{E') , (18) 

essentially the same as in Ref. [l(| at low energy, how- 
ever, this form has a smoother transition to the high 
energy behavior. 

We calculate the Z-moments for charm production tak- 
ing into account the effects of parton saturation at high 
energies as in Ref. 26]. These moments are energy de- 
pendent, with weak dependence at higher energies. In 
our calculation of the neutrino flux, for simplicity we use 
the constant value of the Z-moments for charm given in 
Table [TTJ We neglect the production of other charmed 
particles. We estimate that the inclusion of A c produc- 
tion would increase the flux of neutrinos from charm by 
~ 15% [I?} . For 7T^~ and K ± production, however, we 
assume Feynman scaling, leading to energy independent 
Z-moments. To compute the Z-moments for pp we use 
the parametrization of the rapidity distribution [28| 



^=0.12^-^" 



(19) 



where xe = E^/Ep, and for pj we fit this form to the 
data for tt° production in Ref. [29| . For A' ± we take the 
7T^~ result rescaled by 0.1 [3(|. The decay Z-moments are 
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FIG. 3: Muon neutrino plus antineutrino flux from a slow jet 
supernova in the observer's frame, without including neutrino 
oscillation effects as they propagate on their way to Earth. 
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FIG. 4: Muon neutrino plus antineutrino flux from a gamma 
ray burst in the observer's frame, without including neutrino 
oscillation effects as they propagate on their way to Earth. 



calculated using the expressions in Refs. [3l|, [32J. The 
values for Z-moments are shown in Table [TTJ 

To evaluate the neutrino flux from a source in the 
Earth's frame, co-moving energies are scaled by the jet 
bulk Lorentz factor Tj. The neutrino flux is related to 
the jet luminosity Lj, the distance to the source (II, taken 
here to be 20 Mpc = 6.17 x 10 25 cm, and the jet opening 
angle shown in Table I. The neutrino flux is 
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reff 



t dec/ fhad 
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had\ 



x (Z NM £l + Zl M £ h N ad ) 
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-E~ 
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TABLE II: Z-moments and decay lengths used in the evalua- 
tion of the muon neutrino plus antineutrino flux. The p — > p 
moment is taken to be Z^ffx = 0.5. 



M 


Zmv 


Znm 


77 


ct [cm] 




0.061 


0.55 


0.13 


780 


K± 


0.19 


0.055 


0.0065 


370 


D ± 


0.045 


2.4 x 10" 3 




3.2 x 10" 2 


D° 


0.017 


5.6 x 10" 3 




1.2 x 10" 2 



celeration cross over, and E' min = m p c 2 is the proton's 
rest energy. In Fig. [3] we show the neutrino flux, <f> v , as 
a function of neutrino energy for the slow jet supernova 
model, and in Fig. 01 for the gamma ray burst model. 
The fluxes shown are from pp and pry interactions. At 
low energies, the dominant sources of neutrinos are pions 
and kaons produced in pp interactions. Above approxi- 
mately E ~ 10 GeV, P7 production of pions and kaons 
becomes comparable in strength for the GRB, while for 
the SJS it is subdominant for all energies. 

At low energies, where cooling is not important for the 
D mesons, the D + and D° contributions are approxi- 
mately equal using the parameters in Table [TTJ At high 
energies, pp production of charmed mesons, with their 
shorter lifetimes, dominate although the maximum pro- 
ton energy in the GRB jet restricts the dominance of the 
neutrino flux from charm to a narrow energy interval. 
As noted in Eq. (JSJ), the acceleration time depends on 
unknowns including k. If k = 1, the maximum proton 
energy in the comoving frame increases by a factor of 10, 
increasing the range of energies for which charm decays 
contribute to the neutrino flux. 

To discuss the potential for neutrino flux attenuation 
in the source, we review the proton and neutrino optical 
depths for the slow jet supernova model, where the jet 
does not emerge from the source. The jet propagates to 
r j ~ 5 x 10 10 cm in the SJS model. The optical depth 
t' = an'rj/Tb for protons in pp and pry interactions are 

r' pp ~ 3 x 10 5 

t' v1 ~ 10 6 above threshold 

The proton-proton cross section is a factor of 10 s larger 
than the neutrino cross section with protons at E v = 10 6 
GeV. At high energies, the cross section scales like E® A4 , 
so even at E„ — 10 8 GeV, the proton cross section is a 
factor of more than 10 7 larger than the neutrino cross 
section. At 10 6 GeV, the optical depth for the neutrino 
within the jet is 

t' vv ~ 2 x 10- 2 . 

Neutrino attenuation within the jet is therefore not im- 
portant. 

The densities inside and outside of the jet are not the 
same. In the SJS model, the jet does not reach the stel- 
lar envelope, however, Tj is large enough that the relevant 



The quantity -Ej nax is the energy where cooling and ac- densities would be for He or H. Taking as an example [3c 
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p = (rj/r) 3 g/cm 3 for the He density, the optical depth 
for a 10 6 GeV neutrino is r vp ~ 10. On the other hand, 
if the relevant density is for the H envelope, the den- 
sity is much lower and the optical depth r up <C 1. This 
would indicate that some attenuation may be important 
in a full evaluation of neutrino fluxes for these types of 
sources with a more detailed model of the stellar exterior 
[l4j|. For GRBs, the hadronic jet emerges, so given that 
the neutrino-nucleon cross section is so much lower than 
the proton-proton cross section, neutrino attenuation in 
GRBs should not be important. 

Taking into account the uncertainty in the acceleration 
differential index does not have a significant effect on the 
relative importance of the various contributions to the 
neutrino flux. The high energy neutrino flux in the two 
examples shown here is dominated by charmed meson 
decay. In contrast to pions and kaons, charmed mesons 
give equal amounts of v e + D e and + at the source 
since the neutrinos are coming from semileptonic decays. 
This could give different flavor ratios of neutrinos in neu- 
trino telescopes at very high energies. This has also been 
noted in Ref . [|[ . Since many properties of astrophysical 



sources are not well understood, neutrino measurements 
may hold the key to understanding the environments in 
which they were produced. Charm production in astro- 
physical jets provides an important enhancement of the 
neutrino flux at very high energies. 

Current neutrino experiments, such as IceCube which 
has only partial implementation of the strings, have al- 
ready put some limits on neutrino flux from a specific 
astrophysical source [34(. Once all the strings are imple- 
mented in IceCube, the sensitivity for detection of the 
astrophysical neutrinos will improve significantly. In ad- 
dition, future cubic kilometer deep-sea neutrino detector 
[ol| will also have a very good chance to look for astro- 
physical neutrinos from point sources. 
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